CONVERGENCE AND STABILITY OF THE INVERSE SCATTERING 

SERIES FOR DIFFUSE WAVES 



SHARI MOSKOW AND JOHN C. SCHOTLAND 

Abstract. We analyze the inverse scattering series for diffuse waves in random media. 
In previous work the inverse series was used to develop fast, direct image reconstruction 
algorithms in optical tomography. Here we characterize the convergence, stability and 
approximation error of the series. 



1. Introduction 

The inverse scattering problem (ISP) for diffuse waves consists of recovering the spatially- 
varying absorption of the interior of a bounded domain from measurements taken on its 
boundary. The problem has been widely studied in the context of optical tomography — an 
emerging biomedical imaging modality which uses near-infrared light as a probe of tissue 
structure and function [1, 17]. More generally, diffusion of multiply-scattered waves is a 
nearly universal feature of wave propagation in random media [19]. 

The ISP in optical tomography is usually formulated as a nonlinear optimization prob- 
lem. At present, the iterative methods which are used to solve this problem are not well 
understood mathematically, since error estimates and convergence results are not known. In 
this paper we will show that, to some extent, it is possible to fill this gap. In particular, we 
will study the solution to the ISP which arises from inversion of the Born series. In previous 
work we have utilized such expansions as tools to develop fast, direct image reconstruction 
algorithms [10]. Here we characterize their convergence, stability and reconstruction error. 

Early work on series solutions to the quantum-mechanical inverse backscattering problem 
was carried out by Jost and Kohn [7], Moses [13] and Prosser [16]. There has also been 
more recent work on the inverse medium problem for acoustic waves [5, 18,20]. However, 
the procedures employed in these papers are purely formal [14]. 

The remainder of this paper is organized as follows. In Section 2 we develop the scattering 
theory of diffuse waves in an inhomogeneous medium — this corresponds to the forward 
problem of optical tomography. We then derive various estimates which are later used to 
study the convergence of the Born series and its inverse. The inversion of the Born series is 
taken up in Section 3, where we also obtain our main results on the convergence, stability 
and approximation error of the inverse series. In Section 4 we discuss numerical applications 
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of our results. Finally, in Section 5 we show that our results extend to the case of the ISP 
for propagating scalar waves. 

2. Forward Problem 

Let J] be a bounded domain in R n , for n > 2, with a smooth boundary dil. We consider 
the propagation of a diffuse wave in an absorbing medium whose energy density u satisfies 
the time-independent diffusion equation 

- V 2 u(x) + k 2 {\ + j]{x))u(x) = , ieO, (2.1) 

where the diffuse wave number k is a positive constant and rj(x) > —1 for all x € f2. The 
function r/ is the spatially varying part of the absorption coefficient which is assumed to 
be supported in a closed ball B a of radius a. The energy density is also taken to obey the 
boundary condition 

u{x) + Iv ■ V«(s) = , x £ dfl , (2.2) 

where v is the unit outward normal to d£l and the extrapolation length £ is a nonnegative 
constant. Note that k and rj are related to the absorption and reduced scattering coefficients 
[i a and y! 8 by k = s/3fi^7 s and r](x) = 5fi a (x)/p, a , where 6(i a (x) = p, a - Ma 0*0 and fl a is 
constant [12]. 

The forward problem of optical tomography is to determine the energy density u for a 
given absorption rj. If the medium is illuminated by a point source, the solution to the 
forward problem is given by the integral equation 

u {l) ^)-^fG MH y )d y, X eU. (2.3) 
Jo, 

Here itj is the energy density of the incident diffuse wave which obeys the equation 

- V 2 Ui(x) + k 2 Ui(x) = 6(x - xi) , x € Q , x\ G 30 (2.4) 

and G is the Green's function for the operator —V 2 + k 2 , where G obeys the boundary 
condition (12. 2\ . We may express G as the sum of the fundamental solution 

e -k\x-y\ 

G {x, y) = — , x,yen (2.5) 

4:7r\x — y\ 

and a solution to the boundary value problem 

- V 2 F(x) + k 2 F(x) = , i£fl (2.6) 
F(x)+fi/-VF(x) = -G (x,y) - £v ■ VG (x,y) , x e dn , 

for each |/6 0. That is, we have G = Go + F. By [6] Theorem 6.31, -F € C 2 (0) when y is 
in the interior of 0. 
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The integral equation (12. 3D has a unique solution. If we apply fixed point iteration 
(beginning with itj), we obtain an infinite series for u of the form 



u(x) = Ui{x) — k / G(x,y)r/(y)ui(y)dy 
Jn 

+ k A [ G(x, y)r ] (y)G(y, y') V (y') Ul (y')dydy' + ■■■ . (2.7) 



We will refer to (|2.7p as the Born series and the approximation to u that results from 
retaining only the linear term in 77 as the Born approximation. 

It will prove useful to express the Born series as a formal power series in tensor powers 
of 77 of the form 

4> = Kir] + K 2 7] ®ri + K 3 r] <gi 77 (8> 77 H , (2.8) 

where <j) = U{ — u. Physically, the scattering data (f)(xi, x^) is proportional to the change in 
intensity measured by a point detector at x 2 € dil due to a point source at x\ £ <9f2 [10]. 
Each term in the series is multilinear in 77 and the operator Kj is defined by 

(Kjf) (x 1 ,x 2 ) 

= (-\y +1 k 2:i / G(xi,y 1 )G(yi,y 2 ) ■ ■ ■ G(y j ^ 1 ,yj)G(y j ,x 2 )f{yi, ■ . . ,yj)dyt ■ ■ ■ dyj , 

JBaX — XBa 

(2.9) 

for x\,x 2 € dCl. To study the convergence of the series (12. 8h we require an estimate on the 
L°° norm of the operator Kj. Let / G L°°(B a x • • • x B a ). Then 

H^i/llL°°(9nx9C) = sup |(^j/)(^l,^2)| (2.10) 

< ll/lloo sup k 2j I \G(x 1 ,y 1 )---G(y j ,x 2 )\dy 1 ---dy j . 

We begin by estimating the above integral for j = 1, 

||-Kl||oQ < sup k 2 I \G(xi,y)G(y,x 2 )\dy 

(x 1 ,x 2 )edQxdn JB a 

< k 2 \B a \ sup sup \G(x,y)\ 2 . (2.11) 

x£B a y£dU 

Now, for j > 2, we take out the first and last factors of G in the integral to obtain 
\\Kj\\oo < sup sup \G(xi,y 1 )G(y j ,x 2 }\ 

(xi,X2)e9nx90 yi&B a ,yj£B a 



k 2j [ \G(y 1 ,y 2 )---G(y j - 1 ,y j )\dy 1 ---dy j . (2.12) 

J B„x---xB n 



We then rewrite this as 



\Kj\\oo < sup sup \G(x,y)\) Ij-i , (2.13) 

v x&B a y&dn 



We also note that 
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where 

Ij-i = k 2j [ \G( yi ,y 2 ) ■ ■ • y 3 )\d yi ■ ■ ■ d Vj . (2.14) 

JB«x-xB„ 

Next, we estimate Ij-i recursively: 
Ij-i < sup k 2 \G(yj-x,yj)\dyj 

yj-l&B a J B a 

■k 2 ^ 2 [ \G(y 1 ,y 2 )---G(y j - 2 ,y j - 1 )\dy 1 ---dy j _ 1 , (2.15) 

JB a X — xB a 

from which it follows that 

Ij-l < HooIj-2 , (2.16) 

where 

Poo = sup k 2 \\G(x,-)\\ L i {Ba) . (2.17) 

x€B a 

h = k 4 I' \G(x,y)\dxdy (2.18) 

JBaXBa 

"Si k \B a \\i 00 . 

Ij-i < k 2 \B a \^ . (2.19) 
z/oo = k 2 \B a \ sup sup \G(x,y)\ 2 . (2.20) 

x€B a y£dQ 

Note that the smoothness of solutions to (j2.6f) implies that /Xoo and ^oo are bounded. Making 
use of (|2.12[) and (|2.19p we obtain the following lemma. 

Lemma 2.1. The operator 

Kj : L°°(B a X • • • X B a ) — ► L°°(dn x dtt) 

defined by \2. 9\) is bounded and 

where /Xqq and are defined by j2.17\ ) and \2. 20\) . respectively. 

We now obtain similar L 2 estimates on the norm of Kj. Let / G L 2 (B a x • • • x B a ). Then 



Thus 
Define 



ll*i/llz»(«lx«l) = / |(^/)(x 1 ,x 2 )|^x 1 dx 2 . (2.21) 
From (|2.9p and the Cauchy-Schwarz inequality we have 

\(K j f)(x 1 ,x 2 )\ < k 2 i\\f\\ L 2 ( [ \G(x 1 ,y 1 )G(y 1 ,y 2 )---G(y v x 2 )\ 2 dy 1 .--dyX /2 . 

y JB a X-XBa ' 



Convergence and Stability of the Inverse Scattering Series for Diffuse Waves 

We begin by estimating ||-Ki||: 



k 2 ( / / \G(xi,y)G(y,x 2 )\ 2 dx 1 dx 2 dy 

\Janxdu JB a 



1/2 



< k 2 \B a \ 1/2 sup ||G(x,-)||£ 2(an) . (2.22) 



Now, for j > 2, we find that 



A\ll> sup / \G(xi,yi)G(x2,yj)\ 2 dxidx2 

B a ,yj&B a JdQxdfl 



L jll2 



■ [ \G{yuy 2 ) ■ ■ ■ G( yj ^ y] )\ 2 d yi ■ ■ ■ d Vj , (2.23) 

JB a x — xB a 

which yields 

\\K,\\ 2 < ( sup \\G{x,-)\\ L 2 {dn) ) 2 J 2 ^ , (2.24) 

K x&B a ' 

where 

J i-i = k4J [ \G(yi,V2) ■ ■ ■ G{y j - U y j )\ 2 dy 1 • ■ • dy 3 . (2.25) 

JB a x — xB a 

We now show that Jj-i is bounded, which implies that the kernel of Kj is in L 2 (dft x d£l x 

B a x • • • x B a ) and thus Kj is a compact operator. To proceed, we estimate Jj-i as follows. 

J|_!< sup fc 4 / |G( % --i )% -)| 2 ^ fc4i ~ 4 

y j - 1 £B a JB a 

■I \G(y 1 ,y 2 )---G(y^ 2 ,y^ 1 )\ 2 dy 1 ---dy j - 1 . (2.26) 

Thus, we find that 

Jj-l < Jj-2 , (2-27) 

where 

M2 = sup fc 2 ||G(x,-)||L2(B o) • (2.28) 

Noting that J\ < k 2 \B a \ 1 / 2 ijl 2 , we obtain 

Jj_i < fc 2 !^! 1 / 2 ^ -1 . (2.29) 

Define 

^ 2 = fc 2 |i?a| 1/2 sup ||G(x, .)!&(«,) • ( 2 - 3 °) 

Note that by the smoothness of solutions to (|2.6p and v 2 can be seen to be bounded. 
Using (|2.24p we have shown the following. 
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Lemma 2.2. The operator 

Kj : L 2 (B a x • • • x B a ) — ► L 2 (dQ x dQ) 
defined by \2. 9\) is bounded and 

\\Kj\\2 < V2V?2~ l ' 
where \ii and v-i are defined by 112. 28\) and 112. 30\) . respectively. 

It is possible to interpolate between I? and L°° by making use of the Riesz-Thorin 
theorem [8] . If < a < 1 then 

\\KjW 2 < WKjWI^WKjW^ . (2.31) 

1 — a 

From Lemmas 12.11 and 12.21 we have that 



2 



= {pflZtf 1 W""^) ■ (2-32) 
We thus obtain the following lemma. 
Lemma 2.3. The operator 

Kj : L p (B a x • • • x B a ) — ► L p (dQ x dQ) 
defined by \2. 9\) is bounded for 2 < p < oo and 

lip — ^pf^p i 

where for 2 < p < oo 

2 ^_ 2 

and 

2 1—2 

#ere /U2,f 2,^00 «^ are defined by \2. 28\) . 112. 30\) . \2.17\) and (2.20\) . respectively. 

To show that the Born series converges in the LP norm for any 2 < p < 00, we majorize 
the sum 

J2\\K j r)®---®r ] \\ LPimxdn) (2.33) 

j 

by a geometric series: 

^2\\Kjr]0 - - ® r}\\u>(d!lxdCi) < ^2\\ K j\\p\\ 

j j 

< ^E(/*NliW , (2-34) 
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which converges if MplMUp(B a ) < 1- When the Born series converges, we may estimate the 
remainder as follows: 



N 

y] KjT] ® • • • ® 7] 

3=1 



LP(dflxdfl) 



oo 



< 



j=N+l 

oo 

£ 

i=iV+l 



,i\f+l 



/i p 1 — At p ||r/|| iP ( Ba ) 



(2.35) 



We summarize the above as 

Proposition 2.1. If the smallness condition \\rj\\L P m a ) < l/Vp holds, then the Born series 
\2. 8\) converges in the L p {dVL x d£l) norm for 2 < p < oo and the estimate \2. 35}) holds. 



We note that L°° convergence of the Born series for diffuse waves has also been considered 
in [11]. Corresponding results for the L°° convergence of the Born series for acoustic waves 
have been described by Colton and Kress [4]. 



3. Inverse Scattering Series 

The inverse scattering problem is to determine the absorption coefficient n everywhere 
within f2 from measurements of the scattering data <fi on d£l. Towards this end, we make 
the ansatz that rj may be expressed as a series in tensor powers of 4> of the form 

7? = KL x <f> + K 2 <t> ® 4> + K.3<t> ® 4> ® + • , (3.1) 

where the IC^s are operators which are to be determined. To proceed, we substitute the 
expression (|2.8p for eft into (|3.ip and equate terms with the same tensor power of rj. We 
thus obtain the relations 

KxK x = I , (3.2) 
KvKx ®Kx+ K X K 2 = , (3.3) 
ICsKi ®K 1 ®K 1 + K 2 K X ®K 2 + JC 2 K 2 ®K X + K X K Z = , (3.4) 

3-1 

>Cm Yl K^®---® K im + KjK x ® • • • <g> K x = , (3.5) 

m=l iiH Hm=j 
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T] = • + 




+ 



Figure 1. Diagrammatic representation of the inverse scattering series. 



which may be solved for the ICj's with the result 

£1 = K+ , (3.6) 

/C 2 = -K x K 2 Ki (8) /Ci , (3.7) 

£ 3 = -(/C 2 Ki®K 2 + /C 2 iv:2®^i+/Ci^3)/Ci®/Ci®/Ci , (3.8) 

= - ( ^ /C m ^ ifi! ® • • • <8> if im J Kx ® • • • (8) £i . (3.9) 

\m=l iiH him=j / 



We will refer to (13. ip as the inverse scattering series. Here we note several of its properties. 
First, Kt is the regularized pseudoinverse of the operator K%, not a true inverse, so (|3.2p is 
not satisfied exactly. The singular value decomposition of the operator can be computed 
analytically for special geometries [12]. Since the operator K\ is unbounded, regularization 
of K x is required to control the ill-posedness of the inverse problem. Second, the coefficients 
in the inverse series have a recursive structure. The operator ICj is determined by the 
coefficients of the Born series K\, K2, ■ ■ ■ ,Kj. Third, the inverse scattering series can be 
represented in diagrammatic form as shown in Figure [TJ A solid line corresponds to a factor 
of G, a wavy line to the incident field, a solid vertex (•) to tC\(p, and the effect of the final 
application of K,\ is to join the ends of the diagrams. Note that the recursive structure of 
the series is evident in the diagrammatic expansion which is shown to third order. 



Remark 3.1. Inversion of only the linear term in the Born series is required to compute the 
inverse series to all orders. Thus an ill-posed nonlinear inverse problem is reduced to an ill- 
posed linear inverse problem plus a well-posed nonlinear problem, namely the computation 
of the higher order terms in the series. 



Convergence and Stability of the Inverse Scattering Series for Diffuse Waves 9 

We now proceed to study the convergence of the inverse series. We begin by developing 
the necessary estimates on the norm of the operator ICj. Let 2 < p < oo. Then 

3-1 

ll^jllp — X/ II^Trellpll-^ii lip ' ' ' \\Kim llpll^l lip 

m=l Hm=J 

3-1 

^ H^ilIpE E II^Hp"^ 1-1 ■ • ■ vpifr- 1 , (3-10) 

m=l ixH Hm=j 

where we have used Lemma 12.31 to obtain the second inequality. Next, we define II(j, m) to 
be the number of ordered partitions of the integer j into m parts. It can be seen that 

U(j,m) = ( J ~\) , (3.11) 

3-i 

£ll0\m) = V-i-i . (3.12) 

m=l 

Thus in the diagrammatic representation of the inverse series shown in Figure [H there are 
2 J ~ 1 — 1 diagrams of order j and II(j, m) topologically distinct diagrams of order j and 
degree m with m = 1, . . . , j — 1. It follows that 

3-1 

\%A\p < \\£i\\iJ2 \\IC m \\ P U(j,m)^^- m 

m=l 

< \\jci\& (e ii^Hp) f E mmK4'A 

\m=l / \m=l / 

< fx; \\ic m \\ p ) ( i " 1 )« _1 " m ) 

\m=l / \m=0 ^ ' J 

3-1 

= M^iiipo^p + ^p) jl E n^ m iip • ( 3 - 13 ) 

m=l 

Thus ||/Cj||p is a bounded operator and 

3-1 

\\Kj\\ P < C^ + ^'IIMpE II^Hp • ( 3 - 14 ) 

m=l 

The above estimate for ||/Cj|| p has a recursive structure. It can be seen that 

II^'Hp = CAilh + u p) \\K-i HpFll^i lip > (3.15) 
where, for j > 2, Cj obeys the recursion relation 

C j+1 = Cj + [(/ip + ||/Ci llpPQ , C 2 = 1 . (3.16) 
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It can be seen that 

3-1 

Cj= n (i+[(i*p+»p)\\£i\\pn ■ ( 3 - 17 ) 

m=2 

Evidently, Cj is bounded for all j since 

3-1 

]nCj < ^ln(l + [(fi p + v p )\\JCi\\ p ] m ) 

m=l 
3-1 

m=l 

where the final inequality follows if (fj, p + f p )||/Ci|| p < 1. A more refined calculation using 
the Euler-Maclaurin sum formula gives 

Li 2 (-(/i p + f p ) ll^i Hp) . 1, // , Mir- ii \ / Q1f y> 

lnCj ~ in^ + ^ii^y + 2 ln(( ^ + ^H^IW ' (3 - 19) 

where Li 2 is the dilogarithm function. We have shown the following. 
Lemma 3.1. Let + f»)||/Ci|L < 1 anc? 2 < p < oo . TTien the operator 

iq : L p (dn x x an) — ► l"^) 

defined by \3. 9\) is bounded and 

||p < C(/i p + v p ) j \\K.i \\ p , 
where C = C(/x p ,u p , \\Ki\\ p ) is independent of j. 



Note that (|3.18|) and (|3.19p give explicit bounds for C. 
We are now ready to state our main results. 

Theorem 3.1 (Convergence of the inverse scattering series). The inverse scattering series 
converges in the L v norm for 2 < p < oo if \\JCi\\ p < l/(/^p + v p ) and ||^i<^||lp(b„) < 
l/(/ip + Vp)- Furthermore, the following estimate for the series limit fj holds 

N (/,.,.. mik II II Jill \ N +! 

fj — ^/Cj0 ® • • • <g> <j> 



3=1 



L p (B a ) 1 - {Hp + t / p)||^l||p||0||LP(9nx9Q) 

where C = C{n p ,v p) ||/Ci|| p ) does not depend on N nor on the scattering data (p. 
Proof. The series Y2j K-j'P ® ' ' ' ® converges in norm if 

^||^®---®0|U nBo) < Ell^llfll^Hionxan) ( 3 - 20 ) 



< C^[( M p + ^)||/C 1 ||p||</»|| iP( - " 



'(af7xan)J 
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converges, where the last bound follows from Lemma 13.11 Clearly, the right hand side of 
(|3.20p converges when 

(lip + * / p)||£iUMUp(9nn) < 1 • (3.21) 
To estimate the remainder we consider 

N 



i=l V a> j=N+l 



< C Y, [(M P + ^)l|/Ci||p||^|| iP(9 nxan)] 5 ' , (3.22) 

j=7V+l 

from which the desired result follows. □ 

We now consider the stability of the limit of the inverse scattering series under pertur- 
bations in the scattering data. 

Theorem 3.2 (Stability). Let \\KLi\\ p < l/([i p + v p ) and let <p\ and 02 be scattering data 
for which M\\JC\ \\ p < l/([i p + v p ), where M = max(||0i||p, ||02||p) and2<p< oo. Let rji 
and r\2 denote the corresponding limits of the inverse scattering series. Then the following 
estimate holds 

\\m - m\\LP(B a ) < c\\4>i - 4>2\\LP(dnxdn) , 

where C = C([i p ,u p , ||/Ci|| p ,M) is a constant which is otherwise independent of (pi and fa- 
Proof. We begin with the estimate 

||*7l - m\\LP(B a ) < Y H^i^ 1 ® • • • ® <^l - 02 ® • • • ® fa)\\LP{B a ) ■ (3-23) 
j 

Next, we make use of the identity 

01 ® • • • 01 - 02 ® ■ ■ • ® 02 = 4> 02 • • • <8> 02 + 01 ® 1p 02 ' ' ' 02 

H h 0i ® 0i ® • • • ® V 02 + 01 ® <Pi ® • • • ® 01 ® V> , (3-24) 

where ip = 0i — 02- It follows that 

||?7l - *72 || Z,P(B a ) < ||/Cj|| p ||0l • • • 01 -0 02 ® • • • 02||LP 

j k=l 

= Y^\\p MJ ^\\LP(dnxdn) , (3-25) 
j 

where ip is in the kth position of the tensor product. Using Lemma 13. 11 we have 

Yii(Vp + v P )\\£i\\pM} j 
j 

c 

< Il^l|lpll01 - 02 1 1 L? (Mx Of})- " j 2 ' ( 3 - 26 ) 

[1 - ([i p + fp)||/Ci|| p MJ 

The above series converges when ([i p + f p )||/Ci \\ P M < 1, which holds by hypothesis. □ 



12 Shari Moskow and John C. Schotland 

Remark 3.2. It follows from the proof of Theorem 13.21 that C is proportional to ||/Ci|| p . 
Since regularization sets the scale of ||/Ci|| p , it follows that the stability of the nonlinear 
inverse problem is controlled by the stability of the linear inverse problem. 

The limit of the inverse scattering series does not, in general, coincide with r\. We 
characterize the approximation error as follows. 

Theorem 3.3 (Error characterization). Suppose that ||/Ci|| p < 1 / \(jl p + u p ) , ||^i0[|li»(b o ) < 
+ Up) and 2 < p < oo. Let M. = max (||r/|| iP ( Ba ), \\K.\Kiri\\ LP ^ Ba ^) and assume that 
A4 < l/(/J p + v v ). Then the norm of the difference between the partial sum of the inverse 
series and the true absorption obeys the estimate 



rj- 

3=1 



J2 W ® • • • ® <t> Tpm , < C\\(I - K&Mvm + V 



LHBu)- " V ^n^ a) ~ 1- (Atp + J/^ll/CiHpll^Hp ' 

where C = C(/x p ,u p , \\Ki\\ p ,ftA) and C = C(fj, p , v p \\JCi\\ p ) are independent of N and cj>. 
Proof. The hypotheses imply that the regularized inverse scattering series 

fi = ^K,j(t)®---®(j) (3.27) 

3 

converges. The forward series 

4> = K JV ® • • • ® V (3-28) 

3 

also converges by hypothesis, so we can substitute it into (|3.27p to obtain 

fj = KjV ® • • • ® V , ( 3 - 29 ) 



3 



where 



and 



K\ = ICxKx , (3.30) 
^3 =\J2 }Cm Yl K h ®---®K im j +K j K l ®---®K l , (3.31) 

y m =l ixH h«m=j / 

for j > 2. From (pT9j) it follows that 
j'-i 

^•=^/C m ^ l( g)---®^ m (/-/CiE:i^---(g)/CiKi), (3.32) 

m=l iiH Mm=j 

and so we have 

77 = K\K\r\ + /C 2 r/ ® 77 -\ . (3.33) 

We thus obtain 

r]-rj={I- K x K x )i] - £1^2 (v ® »7 - £i#i»7 ® £i#i?7) + • ■ ■ , (3.34) 
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which yields the estimate 
3-1 

\\n-v\\v -EE E ll^mllpll^illp " " " \\K im \\ p \\ri®- ■ ■®7]-lCiK 1 r]®- ■ -®lCiKi-r]\\ p . 
j m=l h-\ \-im=j 

(3.35) 

Next, we put 

ifi = T) — K,\Kirj (3.36) 
and make use of the identity (|3,24p to obtain 

\\r}®---®r)-1C 1 K 1 r 1 ®---®1C 1 K l r]\\ p <3Mi- 1 \\i)\\ p . (3.37) 

We then have 



j'-l 



\\v -vhp(B a ) < EE E \\ ,c m\\p\\Ki 1 \\p---\\K im \\ p jM J - 

j m=liiH Hm=3 

3-1 

< EE^ i_1 ii^iW' m )<^" ro Wi». ( 3 - 38 ) 

j m=l 

where we have used Lemma 12.31 and H(j,m) denotes the number of ordered partitions of j 
into m parts. Making use of (|3.1ip . we have 



\V-V\\p < ^T / My Mj ' 1 {j2\\ }C rn\\ P ) fEf^ 1 ) 



j,m. j-l-m 

iSp ^p 



j \m=l / \m=0 

3-1 

^ iiv'IIpEE^^^p + ^ii^iIp • ( 3 - 39 ) 

j m=l 

We now apply Lemma 13.11 to obtain 

3-1 

\\V ~ v\\p < CWpE E jM^fjip + v p ) m+i \\lCx\\p , ( 3 -40) 

j m=l 

since the constant C from the lemma is independent of j. Performing the sum over m, we 
have 

ii, - m, < cm, E^-fa + *y> ^Xm-I • (3 - 41) 

which is bounded since A4(fJ, p + u p ) < 1 and (fj, p + f p )||/Ci|| p < 1 by hypothesis. Eq. (|3.4ip 
thus becomes 

||r?-r?|| p <C||(/-/C 1 K 1 ) ?? || p , (3.42) 

where C is a new constant which depends on fj, p ,v p ,M and ||/Ci|| p . Finally, using the 
triangle inequality, we can account for the error which arises from cutting off the remainder 
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of the series. We thus obtain 

TV 



< 



?7 — ^ K.jTJ (g> • • • ® 7] 



LP(B a ) 



LP(Ba) 



\LP(B a ) 



j=N+l 



N+l 



<r r\\(j r Y»ll ^ (Ofg + ^)||^]J|piiriir 

< c\\(i - KiKtinWi,™ + c 1 _ (fip + Vp)m \ p \\4>\\ P 



Evidently, in the above proof, we have shown the following. 



(3.43) 
□ 



Corollary 3.1. Suppose the hypotheses of Theorem \3.3\ hold, then the norm of the difference 
between the sum of the inverse series and the absorption r\ can be bounded by 

\\V ~ vhp(B a ) < C\\(I - /CiKi)r/|| LP(Ba) , 

where fj is the limit of the inverse scattering series and C = C(/J, p , u p , \\Ki\\ p , M) is a 
constant which is otherwise independent ofij. 



Remark 3.3. We note that, as expected, the above corollary shows that regularization of /Ci 
creates an error in the reconstruction of n. For a fixed regularization, the relation K,\K\ = I 
holds on a subspace of L p {B a ) which, in practice, is finite dimensional. By regularizing KL\ 
more weakly, the subspace will become larger, eventually approaching all of L p . However, 
in this instance, the estimate in Theorem 13.31 would not hold since ||/Ci|L would become so 
large that the inverse scattering series would not converge. Nevertheless, Theorem 13.31 does 
describe what can be reconstructed exactly, namely those 77 for which K\K\r\ = I. That is, 
if we know apriori that r] belongs to a particular finite-dimensional subspace of L p , we can 
choose K,\ to be a true inverse on this subspace. Then, if ||/Ci|| p and ||/Ci0||lp are sufficiently 
small, the inverse series will recover 77 exactly. 

4. Numerical Results 

It is straightforward to compute the constants jjL p and v p for M 3 . In this case G(x,y) = 
Go(x, y) for x,y £ Q, and the measurements are carried out on dQ where QcM 3 . We then 
have 

£2 p p -k\x\ 

4tt 



-dx 



B a 



1 - (1 + ka)e 



-ka 



The condition for L°° convergence of the Born series becomes 

1 



\V\\L°° < 



1 - (1 + ka)e 



—ka 



(4.1) 
(4.2) 
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Figure 2. Radius of convergence of the Born series in the I? 
and L°° ( ) norms. 



We note that when ka <C 1, we have H^Hl 00 ^ 0(l/(ka) 2 ). In the opposite limit, when 
ka S> 1, we obtain ||^||l°° < 1 and the radius of convergence is asymptotically independent 
of ka. 

For the L 2 case we obtain 

1/2 



f'2 



fe2 
47T 



k 2 e -ka/2 



-2k\x\ x 



sinh(fea) 



Thus the Born series converges in L 2 if 



D ka/2 



ML? < 



4irk 



ink 



1/2 



k 2 \sinh(A;a) 



1/2 



(4.3) 



(4.4) 



When ka <C 1, we have ||^||l2 < 0(l/(ka) 2 ), which is similar to the L°° estimate. However, 
when b > 1 we find that ||^||i2 < 0(l/(fca) 3 / 2 ), which is markedly different than the L°° 
result. The dependence of the radius of convergence on ka is shown in Figure [5J 

We now give an application of Theorem 13.11 Making use of the Green's function (|2.5p . 
we have 



< k 2 \B a \ 



v 2 < k z \dn\\B t 



e -2kdist(dn,B a ) 

(4irdist(dn,B a )) 2 ' 

e -2fcdist(dQ,B a ) 



1 1/2. 



(47rdist(dft,£ a )) 2 



(4.5) 
(4.6) 



16 



Shari Moskow and John C. Schotland 



R 



3.5 f 

3 
2.5 

2 
1.5 

1 

0.5 




\ 



V 



10 



ka 



Figure 3. Radius of convergence of the inverse scattering series in the L 2 
( ) and L°° ( ) norms. 



Note that v p is exponentially small. We have computed the radius of convergence R p = 
l/O-tp + i/p) °f the inverse scattering series when is a ball of radius 2a, where a is the radius 
of the ball containing the support of r\ and both balls are concentric. The dependence of 
R p on ka is shown in Figure EJ It can be seen that < 1 and R 2 < 0(l/(fca) 3 / 2 ) when 
ka > 1. 



5. Scalar Waves 

In this section we will analyze the inverse scattering series for the case of scalar waves. 
We will see that although the algebraic structure of the inverse Born series for diffuse 
waves is similar to that for propagating scalar waves, its analytic structure is different. 
This distinction is a reflection of the underlying physical difference between the short-range 
propagation of diffuse waves and the long-range propagation of scalar waves. 

Consider the time-independent wave equation 

V 2 u(x) + k 2 (l + 7](x))u(x) =0 , xeR 3 (5.1) 

for the propagation of a scalar wave field u. The wave number k is a positive constant and 
rj(x) is assumed to be supported in a closed ball B a of radius a, with rj(x) > —1. The field 
u is taken to obey the outgoing Sommerfeld radiation condition. The solution to (15. ID is 
given by the integral equation 



u(x) = Ui{x) + k 2 G(x, y)u{y)ri(y)dy , 

JR3 



(5.2) 
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Figure 4. Radius of convergence of the inverse scattering series in the L°° 
norm for diffuse ( ) and propagating waves ( ). 

where U{ is the incident field and the Green's function is given by 

„ik\x— y\ 

G(x,y) = — - (5.3) 

4:Tr\x — y\ 

The convergence, stability and approximation error of the inverse scattering series corre- 
sponding to (|5.ip is also governed by Theorems 13. 11 13T21 and [3T31 However, the values of the 
constants [i v and v v must be modified. The modifications are a reflection of the physical 
difference between the oscillatory nature of the Green's function (|5.3p and the exponentially 
decaying diffusion Green's function. It can be seen that 



k 2 \B n 



(5.4) 

~ - (47Tdist(dft,5 a )) 2 ' (5 ' 5) 

Thus the radius of convergence is = 0(1/ (ka) 2 ) for ka ^> 1. This means that for small 
scatterers, the inverse scattering series has a large radius of convergence, in contrast to the 
case of diffuse waves where Roo = 0(1) even for large scatterers. The dependence of R^o 
on ka is illustrated in Figure 4 for both diffuse and propagating scalar waves. 

6. Discussion 

A few final remarks are necessary. Numerical evidence suggests that the estimates we 
have obtained on the convergence of the inverse series are conservative since ||/Ci|| p ^> 1 in 
practice [10,12]. Nevertheless, insight into the structure of the ISP has been obtained. In 
particular, the inverse series is well suited to the study of waves which do not propagate over 
large scales such as diffuse waves in random media and evanescent electromagnetic waves 
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in nanoscale systems [2,3]. In the latter case, the inverse scattering series has also been 
shown to be computationally effective [15]. To analyze this problem, our methods must be 
extended to treat the Maxwell equations. Other areas of future interest include the study 
of the inverse Bremmer series [9]. 
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